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ABSTRACT 

AG Dra is a symbiotic variable consisting of a metal poor, yellow giant mass 
donor under-filling its Roche lobe, and a hot accreting white dwarf, possibly sur- 
rounded by an optically thick, bright accretion disk which could be present from 
wind accretion. We constructed NLTE synthetic spectral models for white dwarf 
spectra and optically thick accretion disk spectra to model a FUSE spectrum 
of AG Dra, obtained when the hot component is viewed in front of the yellow 
giant. The spectrum has been de-reddened (E(B-V) = 0.05) and the model fit- 
ting carried out, with the distance regarded as a free parameter, but required to 
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be larger than the Hipparcos lower limit of 1 kpc. We find that the best-fitting 
model is a bare accreting white dwarf with M W( i = O.4M , T e ff — 80,000K and 
a model-derived distance of 1543 pc. Higher temperatures are ruled out due to 
excess flux at the shortest wavelengths while a lower temperature decreases the 
distance below 1 kpc. Any accretion disk which might be present is a only a mi- 
nor contributor to the FUV flux. This raises the possibility that the soft X-rays 
originate from a very hot boundary layer between a putative accretion disk and 
the accreting star. 

Key words, stars: binaries: symbiotic stars: fundamental parameters X-rays: binaries 
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1. Introduction 

AG Dra is a D-Type Symbiotic with an orbital period P or b ~ 550 days. The mass donor 
is a metal-deficient KO-4 giant which has also been classified as a Barium star. Several 
estimates of the mass loss from the donor giant have appeared ranging from 10~ 7 M & /yr 
to 2.5 x 1(T 7 M Q /yr (cf. Mikolajewska et al. 1995; Tomov et al. 2000; Tomov & Tomov 
2002). A key question is how much of the outflow from the cool giant is accreted onto 
the hot component. The hot component is presumed to be a white dwarf though no direct 
detection has been reported, and it remains unclear whether there is an accretion disk (from 
wind accretion) surrounding the hot component. The outburst recurrence time of AG Dra is 
every ~ 15 years and the duration of the outbursts is ~ 3 — 6 years. Hipparcos has yielded 
a lower limit to the distance of d > 1 kpc. Mikolajewska et al's (1995) distance of 2.5 kpc 
can be compared with distance estimates based upon estimates for the cool giant radius 
from its rotational velocity. In particular, Fekel et al. (2003) derived v rot smi = 3.6 ± 1.0 
km/s. Assuming corotation with the orbital period and using the inclination derived from 
the spectropolarimetry, an estimate of the red giant radius, R g = 0.209(+0.151/-0.075) 
a.u., is obtained, which combined with the near-IR K-magnitude = 6.2, corresponds to 
d = 1.5(±0.5)kpc. The errors are set by the combined errors of the rotational velocity and 
inclination estimates. It also appears (Mikolajewska et al. 1995; Smith et al. 1996) that the 
giant donor star in AG Dra is brighter than typical for K III giant stars of solar metallicity. 

AG Dra is not an eclipsing symbiotic system, therefore the orbital inclination is below 
i = 80° and estimates between i = 70° and i = 40° have appeared in the literature by Miko- 
lajewska et al. (1995) and Schmid & Schild 1997), with the latter authors' estimate based 
on spectropolarimetry of the Raman scattered O VI lines. In our model fitting described in 
section 3, we ruled out i = 60° and adopted i = 41°. The reddening has been estimated to be 
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E(B-V) = 0.05 (Mikolawjewska et al. 1995), however, reddening values as high as E(B-V) = 
0.10 have not been ruled out (Young et al. 1995). AG Dra has been classified as a supersoft 
X-ray source (Greiner et al. 1997). It has been established that the X-ray/FUV flux level 
of AG Dra is anti-correlated with the optical maxima (Greiner et al. 1997). 

The mass of the hot component has been estimated to be low. Mikolajewska et al. 
(1995) found M h = 0.4 to 0.6 M for a distance of 2.5kpc, while Tomov et al. (2000) gives 
0.5 M . The hot component is presumably a white dwarf though Mikolajewska et al. (1995) 
suggests it is a sub-dwarf with R w d 0.06 — O.O8i? . Note that a low mass for AG Dra's 
hot component would be consistent with the finding of Muerset et al.(1991) that the masses 
of the hot components in symbiotic systems tend to be low. 

However, all previous determinations of the temperature of the hot component have 
been carried out primarily with either blackbody fits to the X-ray and IUE ultraviolet data 
or by the application of the modified Zanstra method, using the emission line strengths. An 
Analysis of the X-ray data has yielded ~ 160, 000K ±20,000K from ROSAT data (Greiner et 
al. 1997) while Mikolajewska et al.(1995) found T h = 80,000K to 150,000K from blackbody 
fits to the IUE FUV data. Gonzalez- Riestra et al.(1999) found an average Zanstra temper- 
ature of 109,600K ±5400K. Likewise, estimates of the accretion rate have been made using 
IUE data by Mikolajewska et al.(1995), Gonzalez-Riestra et al. (1999) during quiescence, 
Viotti et al. (1999) used the N V (1240) emission measure and by Tomov et al.(2000). 

The acquisition of FUSE spectra by Young et al.(1995) and X-ray data with XMM- 
Newton (Gonzalez-Riestra et al.2008) has offered deepened insight into the hot component. 
There has never been, to our knowledge, a direct comparison between the FUSE spectrum of 
AG Dra and state of the art models of optically thick accretion disks with vertical structure 
or with NLTE white dwarf photospheric models. Moreover, even if the continuum (and 
absorption line) radiation in AG Dra's FUV spectra does not arise from an accretion disk, 
it is just as important to determine the hot star properties by direct comparison of the FUV 
spectra with NLTE high gravity photospheres. This modeling would offer an important 
comparison with indirect methods (e.g. modified Zanstra method) or black body fits to the 
FUV spectra. Why is this comparison important? Because it remains unclear at this time 
what comprises the hot component in AG Dra. Is it the inner edge of what may be an 
accretion disk or the surface of a bloated accreting white dwarf? Does the accretion of wind 
from a giant star underfilling its Roche lobe lead to an accretion disk plus boundary layer 
formation or not? If a disk is present, then can the spectral energy distribution be fitted by 
a realistic accretion disk model? Furthermore, the rate of accretion onto the hot component, 
whether via a disk or by direct gravitational capture of wind outflow from the donor giant, 
holds the key to understanding the mechanism of AG Dra's two kinds of outbursts, as well 
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as its intervals of low activity and quiescence. 

In section 2, we present the archival FUSE data for AG Dra. 

2. Spectroscopic Observations 

The FUSE spectrum (see Young et al. 1995) of AG Dra near phase 0.5 (S312010200) 
was obtained on 03-16-2000 with a start time of 15:57:00 UT. The spectrum was taken with 
LWRS with an exposure time of 2388 seconds. The identified lines are shown in Fig. 1. 
The FUSE emission lines were discussed at length in Young et al. (1995). When the FUSE 
spectrum was obtained, the white dwarf was being viewed in front of the giant companion 
and during a low brightness state of AG Dra. We de-reddened the spectrum with E(B-V) = 
0.05. 

In preparation for the fitting of the FUSE spectrum at 0.5 phase, we masked all emission 
lines, as well as ISM absorption features (see Fig. 3 for details). 

The volume emission measure of the nebular continuum, n 2 e V/And 2 ~ 5 — 10 x 10 13 cm-5 
from the Balmer H/3 line, for a reasonable electron temperature Te ~ 10000 — 20000K. This 
emission measure is also consistent with the observed mm/submm emission (e.g. Mikola- 
jewska et al. 2003 and and references therein). We have not included a contribution from 
the nebular continuum. By model fitting the spectral range of the FUSE spectrum, any 
contribution from a nebular continuum is negligible (Skopal 2005) since both during the 
quiet phase and active phase, the nebular continuum contributes only longward of 1200A. 

However, we have explored whether the temperature of the hot component is constrained 
by the nebular emission lines, based upon the assumption implicit in the Zanstra method 
that they form from photoionization by the hot component. We ran simulations with the 
code CLOUDY (Ferland 1998) in order to constrain the nature of the source powering the 
nebulosity in AG Dra (B.Sabra 2012, in preparation). We generated simulated emission line 
ratios and compared them to the observed line ratios. We set up a grid of simulated nebulae, 
all with solar abundances, and plotted the fluxes of the resulting emission lines in line ratio 
diagrams. Every nebula in the grid is illuminated by an incident continuum from a black- 
body at a given temperature, and has a given hydrogen density and ionization parameter, 
defined as the ratio of hydrogen ionizing photons to hydrogen density. For AG Dra, we chose 
radiation temperatures of 75,000K, 100,000K, and 125,000K to represent the hot compo- 
nent source in AG Dra. The hydrogen densities values are log nn/cm" 3 = 3,6,9,12. The 
ionization parameter values are \ogU = —4, —3, —2, —1, 0, 1, 2, 3, 4. In Fig.2a,b, we display 
the results of CLOUDY photoionization simulations of AG Dra overplotted with FUSE and 
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IUE observed line ratios as listed in Table 2 of Young et al. (2005), and in Table 4 of Miko- 
lajewska et al. 1995 (JD 4820, phase 1.686 values), respectively. We chose these datasets 
since they contain a large number of emission lines. Fig. 2a (left panel) displays the behavior 
of the low ionization line ratios while Fig. 2b (right panel) displays the high ionization line 
ratios. The observed line ratios are denoted by "X", while squares, rhombii, and crosses, 
and result from incident radiation due to a blackbody at 75,000K, 100,000K, and 125,000K, 
respectively. The change in density has a minimal effect on the simulated line ratios. Its 
effect is even less than that of the blackbody temperature. The sequence in the line ratios 
is primarily due to the changing ionization parameter. The change in the blackbody tem- 
perature does not result in a substantial change in the simulated line ratios, especially those 
involving lower ionization lines. The lower ionization lines show better agreement with the 
observations than the higher ionization lines. The latter could posssibly also be explained 
by shock heating. We find that the emission line strengths and their ratios remain fairly 
constant for given temperatures between 50,000K and 125,000K. Hence, they do not provide 
a sufficiently accurate and strong constraint on the hot star temperature. 

3. White Dwarf Atmosphere and Accretion Disk Model Fitting Results 

We used the NLTE option in TLUSTY (Hubeny 1988) and SYNSPEC (Hubeny and 
Lanz 1995) to construct the NLTE model atmospheres and the grid of optically thick, steady 
state accretion disk models of Wade and Hubeny (1998) for the disk modeling. The incli- 
nation of the system, % = 41°, was kept fixed in the fitting. The mass of the white dwarf is 
assumed to be in the range 0.4 — O.6M (log g — 7 — 8), which is consistent with the lower 
masses inferred for the hot components in other symbiotic systems (Muerset et al.1991). We 
used the mass-radius relation of Wood (1995) from H-rich white dwarf evolutionary mod- 
els. Our use of models with white dwarf surface gravities seems more plausible than taking 
gravities lower than log g = 7, since we expect the accreting star to be fully degenerate 
with an extended envelope rather than having the radius and gravity of a subdwarf. In the 
symbiotic literature, the term subdwarf is often used but is a misnomer. The physically 
correct characterization of the hot component is a high luminosity degenerate star, that is, 
high luminosity white dwarfs with extended atmospheres and H-burning shells (Sion et al. 
1978; Sion & Starrfield 1994). Thus the use of the term subdwarf implies a lower gravity, 
semi-degenerate object and should be avoided. 

The distance, d, from the scale factor (normalization) of a white dwarf model photo- 
sphere fit is given by 

d = 1000 (R wd /R e ) I VS, 
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where S is the scale factor of the fit. The scale factor, normalized to a kiloparsec, can 
be related to the white dwarf radius through: Fx( bs) = 4n(R 2 /d 2 )H\( mode iy The factor An 
comes from the conversion of Eddington flux, H\, resulting from the program SYNSPEC, 
to physical flux. 

For an accretion disk fit, the models are normalized to 100 pc from Earth. Therefore, 
the distance from the scale factor of a model accretion disk fit is 

d = 100 / v 7 ^ 

We require any fitting solution to correspond to a model-derived distance greater than 
1 kpc, the Hipparcos lower limit to AG Dra's distance. Thus, we regard any model-derived 
distance d < 1 kpc to be implausible. Accordingly, we searched for a best fitting solution 
with a low mass white dwarf and hence a larger radius. The increase in the stellar radius 
leads to distances greater than 1 kpc, the Hipparcos lower limit. In order to obtain a distance 
of at least 1000 pc, then the surface temperature cannot be lower than 70,000K. All of the 
models have solar composition. Of course, the giant component has low metallicity rather 
than solar composition but it is unlikely that this composition difference will substantially 
affect our fitting results. 

Subject to the above parameter constraints, we first carried out model accretion disk 
fits to the FUSE spectrum. We found that with M wd = O.4M , M = 5 x 10~ 8 M Q /yr, 
and % = 41°, reasonable parameters for AG Dra, there is a flux deficit in the model at the 
shortest FUSE wavelengths, relative to the observed spectrum. When we tried an inclination 
i = 60°, the scale factor-derived distance decreases below 1000 pc. If we fix the distance 
for the FUSE spectrum of AG Dra with the white dwarf in front of the red giant, with the 
distance of AG Dra fixed at 2.5 kpc, a widely used distance estimate, then an accretion rate 
M = 5 x 10 -8 M /yr and i = 41°, yields a poor fit. 

In view of these results with model accretion disks, we tried fits of single temperature 
NLTE white dwarf photosphere models to the FUSE spectrum. All of the models have solar 
composition. For a 0.4M Q white dwarf fit with T e ff = 70,000K, the resulting distance d = 
992 pc, just below the Hipparcos lower limit. For a O.4M white dwarf but with T e jj = 
100,000K, the model-derived distance is 2 kpc, but there is a model flux excess at the shorter 
wavelengths of the FUSE spectral range relative to the observed spectrum. We found that 
the best-fitting model was a 0.4M Q white dwarf with T e ff = 80,000K ±5000K and a model- 
derived distance of 1543 pc. However, our model photosphere fits to the FUSE spectrum 
for O.5M yielded reasonable fits as well but with a closer distance barely consistent with 
the Hipparcos lower limit. In Fig. 3, we display the three different model fits to the FUSE 
spectrum for a white dwarf mass of O.4M . The dashed line is the accretion disk fit, the 
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dotted line is the white dwarf photosphere fit with T e // = 100,000K and the solid line is the 
white dwarf photosphere fit with T e ff = 80,000K. 

At these high temperatures (>70,000K), a white dwarf photosphere alone provides vir- 
tually all of the FUV radiation from AG Dra. Hydrogen shell burning driven by accretion 
may be providing the luminosity of the hot component and may explain the supersoft X- 
rays. If an accretion disk is present, then it provides only a relatively small fraction of the 
FUV flux while the hot shell burning accreting white dwarf provides > 90% of the FUV 
flux. Since the FUSE spectrum was obtained during a low optical brightness state, it is not 
unexpected that a bare white dwarf provides the best model. 

4. Discussion and Conclusions 

The most surprising result of our analysis is that our best estimate of the surface tem- 
perature of the white dwarf is only 80,000K. These high gravity NLTE photosphere and disk 
model fits are, to our knowledge, the first to be carried out for AG Dra. Our results, from 
the analysis of a FUSE spectrum, strongly suggest that the true photospheric temperature 
of the white dwarf in AG Dra is considerably cooler than the temperatures derived by mod- 
ified Zanstra methods. No increased amount of interstellar or circumbinary reddening could 
change this conclusion and make the white dwarf as hot as indicated by the Zanstra method. 
It is possible that other mechanisms than photoionization of the cool star wind (e.g. collision 
shocks, disk coronae) are contributing to the emission line strengths. The Zanstra method 
also ignores a hot component wind which could render the equivalent widths of He II and H 
lines to be an inappropriate measure of the true hot component temperature. 

Because of the (relatively) low mass of the accreting star and low mass accretion rate 
(it was observed during a low state), it is very unlikely that the WD temperature reaches 
100,000K. However, a hot boundary layer could significantly contribute to the FUV con- 
tinuum if there is a disk. For a mass accretion rate of the order of ~ lO _9 M /yr, the 
boundary layer theory predicts an optically thick boundary layer with a temperature in the 
range ~100,000K (Godon et al. 1995) to ~ 200, 000K (Popham & Narayan 1995). Such a 
boundary layer would produce an FUV flux similar to the one we modeled with a 110,000K 
WD , therefore possibly bringing the distance in the fitting closer to the estimated 2.5kpc. 
At the same time, the boundary layer would also explain the observed soft X-rays. 
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Fig. 1. — Line identifications in the FUSE spectrum of AG Dra, before it has been dered- 
dened. ISM absorption lines have been marked with vertical tick line at mid-height in each 
panel. There are some strong and sharp emission lines of He II, Ne V,VI & VII and O VI. Note 
also the very broad O VI doublet emission feature affecting the hydrogen Ly/3. The Hi lines 
are marked below the upper panel. 
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Fig. 2. — CLOUDY photoionization simulations of AG Dra observed FUV emission line 
ratios. The left panel (Fig. 2a) displays the behavior of the low ionization line ratios while 
the right panel (Fig. 2b) displays the high ionization line ratios. The observed line ratios are 
denoted by 'X'. The simulated line ratios are denoted by red squares, dark blue rhombii, and 
and light blue crosses, and result from incident radiation due to a blackbody at 75,000K, 
100,000K, and 125,000K, respectively. 
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Fig. 3. — We display the three different model fits to the FUSE spectrum. The dashed line is 
the accretion disk fit with M wd = O.4M , M = 5 x 1(T 8 M /yr, and i = 41°, the dotted line 
is the white dwarf photosphere fit with T e ff — 100,000K and the solid line is the white dwarf 
photosphere fit with T e ff — 80,000K. The ISM absorption lines and the sharp emission lines 
have been masked before the fit and are marked in blue. We also masked the broad Ovi 
doublet emission feature. 



